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Abstract. Ultra-Luminous X-ray sources are accreting black holes that might represent strong 
evidence of the Intermediate Mass Black Holes (IMBH) , proposed to exist by theoretical studies 
but with no firm detection (as a class) so far. We analyze the best X-ray timing and spectral 
data from the ULX in NGC 5408 provided by XMM-Newton. The main goal is to study the 
broad-band noise variability of the source. We found an anti-correlation of the fractional root- 
mean square variability versus the intensity of the source, similar to black-hole binaries during 
hard states. 
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1. Introduction 
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+3 . Ultra-Luminous X-ray sources (ULXs) are point like, off-nuclear, extra-galactic sources, 

with total observed X-ray luminosities (Lx^lO 39 ergs" 1 ) higher than the Eddington lu- 
minosity for a stellar-mass black-hole (Lx~10 38 ergs" 1 ). The true nature of these objects 
is still open to debate ( |Feng fc Soria 20Tl] IFender fe Belloni 2012j) as there is still no un- 
ambiguous estimate for the mass of the compact object in these systems. If the emission 
ly-^ is isotropic and the Eddington limit is not violated, ULXs must be fuelled by accretion 
y^Q . onto Intermediate Mass Black Holes (i.e. IMBHs) with masses in the range 100 — 10 5 M Q 



0^ 1 (Colbert fe Mushotzky 1999[). It is possible that ULXs appear very luminous due to a 



combination of moderately high mass, mild beaming and mild super-Eddington emission 
and that ULXs are an inhomogeneous population composed of more than one class. 



ULXs have been studied intensively during the last decades (see Feng & Soria 2011 
for a review). As in the case of Black Hole Binaries (BHBs), some ULXs undergo spec- 
tral transitions from a hard emission-dominated state to a soft emission-dominated state 
(see [Belloni 2010, Bcllo m""2011l for a description of the spectral states in BHBs). Dur- 
ing the hard state the high-energy spectra of ULXs show a power-law spectral shape 
in the 3-8 keV spectral range, together with a high-energy turn-over at 6-7 keV, and 
a soft excess at low energies (e.g. IKaaret eTa l. 2006). This soft excess can be mod- 
elled by emission coming from the accretion disc and is characterized by a low inner 
disc temperature of «0.2keV, as expected if the black holes in these sources are indeed 



IMBHs (jMiller et al. 20031 iMiTler et al. 2004]) . Nevertheless, |Soria (2007)] suggested that 
some ULXs are consistent with black holes accreting at moderate rate with masses of 
«50 - 100M o . 

The study of the timing properties of ULXs represent a promising way to confirm 
the associations/similarities with this class of sources and BHBs. The Power Density 
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Spectra (PDS) of BHBs are usually composed by broad noise components and narrow 
components (quasi periodic oscillations, QPOs). Just a few ULXs show a QPO and their 
detection has been proven to be difficult (jHeil et al. 2009P . In this work we apply a recent 
tool (from Munoz-Darias, Motta & Belloni 2011), the Root Mean Square (rms)-Intensity 
Diagram, which has been proved to be useful to map states in BHBs (without the need 
of any spectral information). 

We studied the currently best available dataset from the ULX in the dwarf galaxy 
NGC 5408, located at a distance of D = 4.8 Mpc (jKarachentsev et al. 20021) . NGC 5408 
X-l is very bright, with a peak X-ray luminosity in the (0.3-10 keV) energy range of 



L x = 2xl0 40 ergs~ 1 . |Strohmayer fc Mushotzky (2009)| found a QPO in its PDS centred 
at 0.01 Hz and inferred a mass for the black hole in the range of 10 3 — 10 4 M . Recently, 



Dheeraj & Strohmayer (2012) studied the timing and spectral properties of NGC 5408 



X-l and have found that the QPO frequency is variable and largely independent on the 
spectral parameters. They suggested that NGC 5408 X-l is accreting in the saturation 
(constancy of the power-law photon index and disc flux with a further increase in the 
QPO centroid frequency) regime a few times observed in BHBs (Vigna rca et al. 2003[ ). 
We performed timing and spectral studies, but focusing on the evolution of the broad- 
band noise rather than on the characteristics of the QPO. 



2. Timing and spectral analysis 

We considered the 6 long (120-130 ks) high-quality observations performed with XMM- 
Newton during 6 years (2006-2011) of NGC 5408 X-l. We limited our analysis to the 
EPIC /pn data, due to its higher time resolution and effective area. We refer to Caballero- 
Garcia et al. (2012, in prep.) for a study of the full (pn+MOS) dataset. 

2.1. Timing 

We performed the analysis of the fast time variability of NGC 5408 X-l in the 1-8 keV 
energy range. The (1-8 keV) pn count rate is 0.2—0.3 ctss^ 1 on average. We excluded from 
the analysis the 8-10 keV energy band because of the lack of signal and we ignored the 0.3- 
1 keV energy range because variability is suppressed in this energy range for NGC 5408 
X-l ([Middleton et al. 201 ip . We produced light curves with a time bin of 0.5 s. We used 
custom software under MATLAB to perform Fast Fourier Transforms (FFT) and produce 
the PDS. A total of 2048 poin ts of the light curve w ere used to perform the FFT. The PDS 
were normalized according to Leahy et al. (1983) and converted to square fractional root 



mean square deviation (rms; BelToni & Hasing er 1990[ ). PDS fitting was carried out with 
the standard XSPEC fitting package by using a one-to-one energy-frequency conversion 
and a unit response. 

We fitted the PDS with a model constituted by a broad Lorentzian centred at zero 
frequency plus a constant (expected to be 2 in the Leahy normalization) to account for 
the Poisson noise. This resulted to be a good fit in all the observations, except for the 
second observation, in which residuals in the form of a QPO are visible around 0.01 Hz. 
The integrated fractional rms in the 0.005-1 Hz band for each observation is reported in 
Tab.[U 

2.2. Spectroscopy 

We produced background-subtracted energy spectra in the 0.3-10 keV energy range, with 
a background taken from an annulus around the target. Because our goal is to study the 
intrinsic emission of the ULX, in this way we expected to minimize the diffuse emission 
from the galaxy. Nevertheless, emission lines are present in the soft X-rays (0.3-1 keV) 



Aperiodic variability of NGC 5408 X-l 3 



Table 1. Results from the timing-spectral study 1 ' 2 
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Power- law 
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41±4 


0.154±0.003 


1.76±0.11 


3.3±0.4 


(3.2±0.1)xl0~ 12 


2 


48±5 


0.156±0.004 


1.60±0.20 


3.1±0.7 


(2.9±0.1)xlCr 12 


3 


35±4 


0.155±0.002 


1.87±0.11 


4.0±0.6 


(3.6±0.1)xl0^ 12 


4 


32±4 


0.157±0.002 


2.00±0.10 


5.2±0.9 


(3.5±0.1)xl0~ 12 


5 


39±3 


0.159±0.002 


1.85±0.12 


3.9±0.7 


(3.3±0.1)xl0~ 12 


6 


41±4 


0.162±0.003 


2.10±0.10 


6.8±1.9 


(3.2±0.1)xl0~ 12 



Note: 

1 The results shown are from the model tbabs Capec+diskpn+cutof f pi) . 

2 Errors arc la. 



energy range, which may indicate the presence of residual diffuse emission from the 
galaxy. We considered the model component apec for the description of this diffuse 
emission plus emission from the inner disc (diskpn model component, with corrections 
for the inner disc emission close to the black hole), using the tbabs absorption model 
component. To fit the high-energy emission, we considered model components allowing 
for high-energy curvature, i.e. either a powerlaw with a high-energy cut-off cutoff pi or 
compTT, the latter with the temperature for the input soft photons tied to the temperature 
of the inner part of the disc. A variable multiplicative constant factor was used to account 
for coordinated variability of the spectra. We fitted the spectra simultaneously, in order 
to tie the parameters related to the diffuse emission (apec model) and the column density 
components, expected to be constant. 

In the model with compTT, the obtained temperature of the electrons in the corona was 
kT = 2 — 3keV and the optical depth r = 4 — 5. In the phenomenological model with 
cutoff pi, the photon indices and the high-energy cut-off found are in the range 1.7 — 2.0 
and 3 — 7 keV, respectively. The two models provided the same statistical description to 
the data and none is statistically favoured. The most relevant results of this analysis are 
shown in the following and in Tab.[TJ 



3. Results and Implications 

We found that the X-ray emission of NGC 5408 X-l is highly variable and that this 
variability changes during the course of the observations (i.e. the fractional rms changes 
from «30% to «50%). This level of variability is comparable to BHBs in the low/hard 
state and to some extreme Seyfert galaxies (Vaugha n et al. 2003[ ). The best model applied 
to the energy-spectra has best-fit parameters that closely resemble those obtained from 
BHBs during the low/hard state (r = 1.7-2.0, kT i , 1 saQ.17keV. IReis et al. 2010j) . except 
for the presence of a low-energy cut-off in the spectra (E = 3 — 7keV), much lower 
than those observed in BHBs (E = 50 — 100 keV). This low value for the cut-off can be 
physically understood as a low value of the temperature of the electrons in the corona and 
has previously reported in the study of a sample of bright ULXs (Gladst one et al. 2009|) . 

We studied possible correlations of the rms with spectral parameters. We found that 
there is an anti-correlation of the total flux with the fractional rms (see details in the cap- 
tion of Fig. [I]), similar to that found in BHBs during the low/hard state flMufioz-Darias, Motta fc Belloni 2011[ ). 
A linear rms-flux correlation in NGC 5408 X-l has already been reported ( Heil fc Vaughan 2010]) , 
but on short time scales (i.e. within single XMM-Newton observations; coinciding with our 
Obs. 1 and 2). Our study shows that there is a fractional rms-flux anti-correlation in the 
longer time-scales, in analogy to what observed in the LHS of BHBs (Munoz-Darias, Motta & Belloni 2011 1. 
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Figure 1. Total unabsorbed flux (in the 0.3-10 keV energy range) versus the fractional rms of 
the variability in the 1-8 keV energy range and 0.005-1 Hz frequency range. 



